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Abstract 

We report measurements of B to pseudoscalar-pseudoscalar decays with at least one rj meson in 
the final state using 140 fb~^ of data collected by the Belle detector at KEKB e^e~ collider. We 
observe the decays of B'^ rjn'^ and B'^ ilK^] the measured branching fractions are B{B'^ — > 
r/7r+) = (4.8t|];^(stat) ± 0.3(sys)) x 10"^ and B{B+ r]K+) = (2.1 ± 0.6(stat) ± 0.2(sys)) x 10"^. 
Their corresponding CP violating asymmetries are measured to be 0.07 it 0.15 it 0.03 for t^tt^ and 
-0.49 ± 0.31 ± 0.07 for rjK^. No significant signals are found for neutral B meson decays. We 
report the following upper limits on branching fractions at the 90% confidence level: B[B^ —>■ 
r]K^) < 2.0 X 10-6, B{B^ r/vrO) < 2.5 x 10"^ and B{B^ r]r}) < 2.0 x lO^^. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 
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Charmless B decays provide a rich sample to understand the B decay dynamics and to 
search for CP violation. An unexpectedly large B — > rj'K branching fraction has stimulated 
much theoretical interest. It was suggested even before the rj'K measurement that two b ^ s 
penguin amplitudes are constructive in S — > rj'K decays but destructive in B —y rjK 
The situation is reversed for B rj'K* and B — > rjK* decays. Experimental results have 
more or less confirmed this picture; however, precise measurements of branching fractions are 
needed to quantitatively understand the contribution of each diagram. It was also pointed 
out that the suppressed penguin amplitudes in the r]K mode may interfere with the CKM 
suppressed b ^ u (tree) amplitude and result in direct CP violation The penguin-tree 
interference may also be large in B^ 'q'n^ j3| and B^ 7771 decays although theoretical 
expectations on the partial rate asymmetry (Acp) could be either positive or negative 
Recently, the BaBar collaboration has observed large negative Aqp values in both rjK^ and 
r]'7i~^, which are ~ 2a away from zero j^. However, more data are needed to verify these 
large CP violating asymmetries. Furthermore, branching fractions and Acp in charmless B 
decays can be used to understand the tree and penguin contributions and provide constraints 
on the third unitarity triangle (p^ (g}. 

In this paper, we report measurements of branching fractions and partial rate asymmetries 
for B ^ rjh decays, where h could be a /C, vr or 77 meson. The partial rate asymmetry is 
defined as: 

_ N{B- ^ f-)-N{B+ ^ n 
N{B~ ^ n + N{B+ ^ 

where N{B~) is the yield for the B^ rjh" decay and N{B^) denotes that of the charge 
conjugate mode. The data sample consists of 152 million BB pairs (140 fb~^) collected 
with the Belle detector at the KEKB e^e" asymmetric-energy (3.5 on 8 GeV) colhder 
operating at the T(45') resonance. 

The Belle detector is a large- solid- angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector (SVD), a 50- layer central drift chamber (CDC), an array of 
aerogel threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scin- 
tillation counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals 
(ECL) located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. 
An iron fiux-return located outside of the coil is instrumented to detect mesons and to 
identify muons (KLM). The detector is described in detail elsewhere 0. 

Candidate neutral pions are selected by requiring the two-photon invariant mass to be 
in the mass window between 115 MeV/c^ and 152 MeV/c^. The momentum of each photon 
is then readjusted, constraining the mass of the photon pair to be the nominal 7r° mass. 
To reduce the low energy photon background, each photon is required to have a minimum 
energy of 50 MeV and the vr^ momentum must be above 250 MeV/c in the laboratory frame. 
Two T] decay channels are considered in this analysis: 77 ^ 77 (7/77) and 7/ Tr+Tr^Tr" (773^). 
In the 7777 reconstruction, each photon is required to have a minimum energy of 50 MeV and 
the energy asymmetry, defined as the energy difference between the two photons divided by 
their energy sum, must be less than 0.9. Furthermore, we remove rj candidates if either one 
of the daughter photons can pair with any other photon to form a 7r° candidate. Candidate 
7/37r mesons are reconstructed by combining a 7r° with a pair of oppositely charged tracks, 
originated from the interaction point (IP). We make the following requirements for the 
invariant mass (in MeV/c^) on the 7/ candidates: 516 < M^^ < 569 MeV/c^ for 7/77 and 
539 < Ms^ < 556 MeV/c^ for 7/3^. An 7/ mass constraint is implemented after the selection 
for each candidate. 
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Charged tracks are required to come from the IP. Charged kaons and pions directly from 
B decays are identified by combining information from the CDC (dE/dx), the TOF and 
the ACC to form a K{tt) likelihood Lk^Lt^). Discrimination between kaons and pions is 
achieved through the likelihood ratio Lk/{Lt^ + Lx)- Charged tracks with likelihood ratios 
greater than 0.6 are regarded as kaons, and less than 0.4 as pions. Furthermore, charged 
tracks that are positively identified as electrons or muons are rejected. K^^ candidates are 
reconstructed from pairs of oppositely charged tracks with invariant mass (M^^) between 
480 to 516 MeV/c^. Each candidate must have a displaced vertex with a flight direction 
consistent with a originating from the IP. 

Candidate B mesons are identified using the beam constrained mass, Mbc = 



-^bcam ^ Pb^ the energy difference, AE = Eb — -Ebcam, where -Ebeam is run-dependent 
and determined from B D^*^Tr events, and Pb and Eb are the momentum and energy of 
the B candidate in the T(4S') rest frame. The resolutions on Mbc and AE are around 3 
MeV/c2 and ~ 20-30 MeV, respectively. Events with Mbc > 5.2 GeV/c^ and \AE\ < 0.3 
GeV are selected for the final analysis. 

The dominant background comes from the e~^e~ — > qq continuum, where q = u,d, s or c. 
To distinguish signal from the jet-like continuum background, event shape variables and the 
B flavor tagging information are employed. We form a Fisher discriminant from seven 
variables that quantify event topology. The Fisher variables include the angle 6t between 
the thrust axis of the B candidate and the thrust axis of the rest of the event, five 



modified Fox- Wolfram moments and a measure of the momentum transverse to the 
event thrust axis {S±) jl^. The probability density functions (PDF) for this discriminant 
and cos Ob, where 6*^ is the angle between the B flight direction and the beam direction in 
the T(4S') rest frame, are obtained using events in the signal Monte Carlo (MC) and data 
with Mbc < 5.26 GeV/c^ for signal and qq background, respectively. These two variables 
are then combined to form a likelihood ratio 71 = Cs/ {Cs + Cqq), where Cs(^gg) is the product 
of signal (qq) probability densities. 

Additional background discrimination is provided by the B flavor tagging. We used 
the standard Belle B tagging package [l^, which gives two outputs: a discrete variable q 
indicating the B flavor and a dilution factor (r) ranging from zero for no flavor information 
and unity for unambiguous flavor assignment. We divide the data into six r regions. The 
continuum suppression is achieved by applying a mode dependent cut on 71 for events in 



each r region based on N^^^/^Ns^^ + Ng^^, where N^^^ is the expected signal from MC 
and Nqq^ denotes the number of background events estimated in data. This 7Z requirement 
retains 58-86% of the signal while reducing 82-96% of the background. From MC all other 
backgrounds are found to be negligible except for the riK~^ <-> rjTr^ reflection, due to K^-7r~^ 
misidentification, and the ?7-ft'*(892)(?7p(770)) feed-down to the r]K{r]7r) modes. We include 
these two components in a fit to extract the signal. 

The signal yields and branching fractions are obtained using an extended unbinned 
maximum- likelihood (ML) fit with input variables Mbc and AE. The likelihood is defined 
as: 

N 

C = exp (- N,) n [ E N,P;{M^o, AE) ], 

j i 

where Nj is the yield of category j (signal, continuum background, reflection, rjK* /rjp) , 
Pj{M\yc, AE) is the probability density for the ith event and A^ is the total number of events. 
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The PDFs of the signal, the reflection background and the rjK* /rjp feed-down are modeled 
with two-dimensional Mbc-A£' smooth functions obtained using MC. The peak positions and 
resolutions in Mbc and AE are adjusted according to the data-MC differences using large 
control samples of i? — > Dtt and D K^tt'ti^ /ti^tt^ decays. The continuum background 
in AE is described by a first or second order polynomial while the Mbc distribution is 
parameterized by an Argus function, /(x) = x\/l — exp [—.^(1 — x^)], where x is Mbc 
divided by half of the total center of mas energy. The continuum PDF is the product of 
an Argus function and a polynomial, where ^ and the coefficients of the polynomial are free 
parameters. Since B — rjK* decays were observed with relatively large branching fractions 
|l4| (~ 20 X 10~^), their feed-down to the rjK modes are fixed from MC in the likelihood 
fit. Since the decay — ^ rjp'^ is experimentally poorly constrained, the amount of this 
background in the tj-k modes is allowed to float in the fit. In the charged B modes, the 
normalizations of the reflections are fixed to expectations based on the B^ — riK~^ and 
B^ — i> T]iT^ branching fractions and ^ vr^ fake rates, measured using 
data. The reflection yield is first estimated with the assumed riK~^ and rjir'^ branching 
fractions and is then recalculated according to our measured branching fractions. No BB 
contributions are considered for the B^ — > 7777 mode. 

In Table |l]for each decay mode we show the measured branching fractions as well as other 
quantities associated with the measurements. The efficiency for each mode is determined 
using MC simulation and corrected for the discrepancy between data and MC using the 
control samples. The only discrpancy we find is the performance of particle identification, 
which results a 4.3% correction for the rjir'^ mode and 1.7% for B rjK^ . The combined 
branching fraction of the two r] decay modes is obtained from a simultaneous likelihood fit 
to all the sub-samples with a common branching fraction. The statistical error in the signal 
yield is taken as the change in the central value when the quantity —2 In C increases by one 
unit from its minimum value. The statistical significance is taken as the square root of the 
difference between the value of — 21n£ for zero signal yield and the minimum value. The 
number of B^B~ and B^^ pairs are assumed to be equal. 

Systematic uncertainties in the fit due to the knowledge of the signal PDFs are estimated 
by performing the fit after varying their peak positions and resolutions by one standard de- 
viation. In the T]K modes, we also vary the expected r]K* feed-down by 1 standard deviation 
to check the yield difference. The quadratic sum of the deviations from the central value 
gives the systematic uncertainty in the fit, which ranges from 3% to 6%. The performance 
of the TZ cut is studied by checking the data-MC efficiency ratio using the B~^ D^tt^ 
control sample. The obtained error is 2.4-3.5%. The systematic errors on the charged track 
reconstruction are estimated to be around 1% using partially reconstructed D* events, and 
verified by comparing the ratio of 77 ^ tt+tt^tt'^ to r/ ^ 77 in data with MC expectations. 
The TT*^ and 77^^ reconstruction efficiency is verified by comparing the vr*^ decay angular dis- 
tribution with the MC prediction, and by measuring the ratio of the branching fractions 
for the two t] decay channels: rj ^ and rj — s> vr^vr^vr^. We assign 3.5% error for the 
7r° and 7/^^ reconstruction. The Kg reconstruction is verified by comparing the efficiency 
ratio of KgTc'^ and K^n^n^. The Kg detection systematic error is 4.4%. 

The uncertainty of number of BB events is 1%. The final systematic error is obtained by 
first summing all correlated errors linearly and then quadratically summing the uncorrelated 
errors. 

Figure n shows the Mbc and AE projections after requiring events to satisfy —0.1 < 
AE < 0.08 GeV (-0.15 < AE < 0.1 GeV for the 7/^^ and 7/7r° modes) and Mbc > 5.27 
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TABLE I: Detection efficiency, product of daughter branching fractions, yield, statistical signifi- 
cance, measured branching fraction, the 90% C.L. upper limit and Acp for the B i]h decays. 
The first errors on columns 4, 6 and 8 are statistic and the second errors are systematic. 



Mode 




Yield 


Sig. 


6(10-6) 


UL(10-6) 


Acp 










8.1 


4.8 ±0.7 ±0.3 




0.07 ±0.15 ±0.03 




23.3 


39.4 


73.4t}|^ ± 2.0 


7.2 


5.3tJ:^ ± 0.3 




0.11 ±0.17 ±0.03 




14.8 


22.6 


19.6t^;? ± 0.7 


4.0 


3.8lJ:^ ±0.3 




r, 11+0.35+0.04 
"•^^-0.33-0.05 










4.0 


2.1 ±0.6 ±0.2 




-0.49 ±0.31 ±0.07 




21.1 


39.4 


28.0t^°(° ± 1.6 


3.5 


2.2l[};^ ± 0.2 




-^Ahtofi ± 0.07 




13.8 


22.6 


7.4t|;^ ± 0.5 


1.8 


1.5tJi±0.2 




n 70+1.03+0.11 

'"-0.76-0.12 










0.4 


0.3;[];?±0.1 


< 2.0 






22.9 


13.5 


-1.9j;|f ±0.3 


- 


-0.4lg:? ± 0.1 








12.2 


7.8 


3.5l|^ ± 0.2 


1.4 


2.4l^| ± ±0.3 














1.9 


1.2 ±0.7 ±0.1 


< 2.5 






17.0 


39.0 


lo 0+8.9+0.8 
-■■"•^-8.0-0.7 


2.5 


1.8l[j| ± 0.2 








11.2 


22.3 


-3.0l||] ± 0.3 




-0.8lJi±0.1 






B^ 7]r] 








1.2 


0.7tL'±0.1 


< 2.0 




^77^77 


16.9 


15.5 


-1.5^?:^ ±0.1 




-0.4+[5;J ± 0.0 






^77%7r 


11.3 


17.8 


7.3;to ± 0-2 


2.3 


2.3;j;^ ±0.2 








7.7 


5.1 


0.3l?:^ ±0.1 


0.2 


0.5l?:^ ±0.1 







GeV/c^, respectively. No significant signals are observed for the neutral B meson modes; 
we set their branching fraction upper limits at the 90% confidence level. The upper limit for 
each mode is determined using the combined likelihood with the reconstruction efficiency 
reduced by its systematic error. We vary the signal PDF and the expected rjK* feed-down 
in the rjK^ mode to compute the likelihood; the largest branching fraction that covers 90% 
of the likelihood area is chosen to be the upper limit. 

Significant signals are observed for charged B decays and we investigate their partial 
rate asymmetries by extracting signal yields separately from the B^ and B~ samples. A 
likelihood fit is performed independently for the two 77 decay modes. The same signal and 
background PDFs used in the branching fraction measurement are applied. The parameters 
of the continuum PDF are fixed according to the branching fraction results. Contributions 
from BB backgrounds are required to be equal for the B^ and B~ samples. Figure El shows 
the Mbc and AE projections. The Acp results for the two rj decay modes are combined 
assuming that the errors are Gaussian. Systematic errors that arise from the knowledge of 
the signal PDF are estimated by varying the peak positions and resolutions. We also check 
the Acp values after varying the amount of the expected rjK* feed-down and the reflection 
background. The BB contributions are allowed to be different for the two samples to obtain 
the systematic error. The largest uncertainty is caused by the reflection. A possible detector 
bias in Acp is studied using B — > Dtt^ decays. The obtained uncertainty is 0.5%. Each 
Acp deviation is added quadratically to provide the total systematic uncertainty. 

In summary, we have observed both B^ rfTi^ and B^ — > rjK^ decays. Their measured 
branching fractions and partial rate asymmetries are summarized in Table EJ We conclude 
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(GeV/c' 



0.25 0.25 

AE (GeV) 



FIG. 1: Mbc and AE projections for (a,b) B+ r?7r+, (c,d) B+ r]K+ , (e,f) B° r]K^ , (g,h) 
rjTT^ and (i,j) B^ — > rjrj decays with the tj^^ and modes combined. Open histograms are 
data, sohd curves are the fit functions, dashed hnes show the continuum contributions and shaded 
histograms are the rjK* /rjp contributions. The small contributions around M^c = 5.28 GeV/c^ 
and AE = ±0.05 GeV in (a)-(d) are the backgrounds from B~^ — > rjTr'^ and B^ rjK'^. 



that the r]7r^ branching fraction is larger than that of r]K^. Our measured — > 7771 
branching fraction is consistent with the BaBar result; however, unlike the large negative 
Acp observed by BaBar, our central value is small and positive but is consistent with no 
asymmetry. For the decay B'^ ^iK^ , our measured branching fraction is 40% lower than 
the BaBar result, corresponding to a 1.3 a deviation. Interestingly, both experiments sug gest 
a large negative Acp value for — > rfK^ ^ which is anticipated by some theories [15[. No 
significant signals are found in the neutral B meson decays and we give their upper limits 
at the 90% confidence level. 

We thank the KEKB group for the excellent operation of the accelerator, the KEK 
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(GeV/c') AE (GeV) 

FIG. 2: Mbc and AE" projections for (a,b) B+ r/7r+, (c,d) B' r/Tr", (e,f) 5+ ^ r?/i:+, and 
(g,h) B~ — > r/X" with the rj^-y and 7^3,^ modes combined. Open histograms are data, solid curves 
are the fit functions, dashed hues show the continuum contributions and shaded histograms are 
the r]K*/r]p contributions. Small curves around Mbc = 5.28 GeV/c^ and AE = ±0.05 GeV are 
the reflection background on B'^ r/vr"*" and B^ 7]K^. 

Science and Technology of India; the BK21 program of the Ministry of Education of Korea 
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